Quantum corrections to the mass of self-dual vortices 
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using heat kernel/generalized zeta function regularization methods. 
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1. In this note we shall compute the one- loop mass 
shift for Abrikosov-Nielsen-Olesen self-dual vortices in 
the Abelian Higgs model. Non-vanishing quantum cor- 
rections to the mass of A^ = 2 supersymmetric vortices 
were reported during the last year in papers [ij and |2]- 
In the second paper, it was found that the central charge 
of the N = 2 SUSY algebra also receives a non- vanishing 
one-loop correction which is exactly equal to the one-loop 
mass shift; thus, one could talk about one-loop BPS satu- 
ration. This latter result fits in a pattern first conjectured 
in |3| and then proved in |j| for supersymmetric kinks. 
Recent work by the authors of the Stony Brook/Viena 
group, [5J, unveils a similar kind of behaviour of super- 
symmetric BPS monopoles in Af = 2 SUSY Yang-Mills 
theory. In this reference, however, it is pointed out that 
(2-|-l)-dimensional SUSY vortices behave not exactly in 
the same way as their (1-1-1)- and (3-|-l)-dimensional 
cousins. One-loop corrections in the vortex case are in 
no way related to an anomaly in the conformal central 
charge, contrarily to the quantum corrections for SUSY 
kinks and monopoles. 

We shall focus, however, on the purely bosonic Abelian 
Higgs model and rely on the heat kernel/generalized zeta 
function regularization method that we developed in ref- 
erences p, H and Q to compute the one- loop shift 
to kink masses. Our approach profits from the high- 
temperature expansion of the heat function, which is 
compatible with Dirichlet boundary conditions in purely 
bosonic theories. In contrast, the application of a sim- 
ilar regularization method to the supersymmetric kink 
requires SUSY friendly boundary conditions, see [21]. We 
shall also encounter more difficulties than in the kink case 
due to the jump from one to two spatial dimensions. 

Defining non-dimensional space-time variables, x^ —^ 
— x^, and fields, (j) —^ vcj) — v{(l)i +i(j>2), A^ -^ vA^, from 
the vacuum expectation value of the Higgs field v and 
the C/(l)-gauge coupling constant e, the action for the 
Abelian Higgs model in (2+l)-dimensions reads: 



with C/(0,(/»*) = !((/)*(/)- 1)^. K = ^ is the only classi- 
cally relevant parameter and measures the ratio between 
the masses of the Higgs and vector particles; A is the 
Higgs field self-coupling. For k, = \ one finds self-dual 
vortices with quantized magnetic flux g = -^, / G Z, 
and mass My = Tr\l\v'^ as the solutions of the first-order 
equations Dicj) ± i-D20 = 0, F12 ± li(j)*(j) - 1) = 0, or, 

(9101 + ^102) T (a2</'2 - A24)l) = 
±(^201 -I- ^202) + (9l02 - Ai(j)i) = 
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with appropriate boundary conditions: 0*0|. 
A0|soo = {di(t) - iAi(j))\s^ = 0, that is, (j)\s^ = e 
and ^ilsoo — ^'i4>*9i<j>\s^- In what follows, we shall fo- 
cus on solutions with positive /: i.e., we shall choose the 
upper signs in the first-order equations. 

2. L^-integrable second-order fluctuations around a 
given vortex solution are still solutions of the first-order 
equations with the same magnetic flux if they belong to 
the kernel of the Dirac-like operator, T>£^{x) — 0, 10] 
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where ^'^(x) ~ {ai{x),a2{x),(pi{x),(p2{x)). We denote 
the vortex solution flelds as ^ = i^i-\-i^2 a-nd V}~, k = 1,2. 
Assembling the small fluctuations around the solution 
4>{x) = 4'{x) + ^{x), Ak{x) = Vk{x) + ak{x) in a four 
column $,{x), the first component of V^ gives the de- 
formation of the vortex equation (3), whereas the third 
and fourth components are due to the respective defor- 
mation of the covariant holomorphy equations (2) and 
(1). The second component sets the background gauge 
Biflk, ^'^Tp) — dkcik — (V'i¥'2 — ip2'fii) on the fluctuations. 
The operators 
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are defined as 7Y+ =15^2? -the second order fluctuation 
operator around the vortex in the background gauge- and 
its partner H^ — VD^ . 

One easily checks that dimkerPt = q Thus, the 
dimension of the moduli space of self-dual vortex solu- 
tions with magnetic charge I is the index of V: indP = 
dimkerl? — dimkerP^. We follow Weinberg |1^], us- 
ing the background instead of the Coulomb gauge, to 
briefly determine indP. The spectra of the operators 
Ti~^ and Ti.~ only differ in the number of eigen-functions 
belonging to their kernels. For topological vortices, we 
do not expect pathologies due to asymmetries between 
the spectral densities of 7i+ and TC~ and thus indP = 
rpj.g-/3w+ _ Tre^'3^". See [ilE3 for the case of Chern- 
Simons-Higgs topological vortices. 

The heat traces Tre~^^ — tr/j^j d^x K-^±{x,x; f3) 
can be obtained from the kernels of the heat equations: 

^I + W±W„±(f,y;/3)=0 

Kn±{x,y;0)^l-S^^\x~y) 

Bearing in mind the structure Ti^ = — A I + I + 
Q^ {x)dk+V^ (x) , one writes the heat kernels in the form: 

Kh± {x, y; 13) = C^{x, y; (3)K-Ho {x, y\ fi) 

with C±(f,x;0) = I. Knoix,y;P) = |^ • I • e"^ 
is the heat kernel for the Klein-Gordon operator Ho = 
(— A-|-1)I, which is the second-order fluctuation operator 
around the vacuum in the Feynman-'t Hooft renormaliz- 
able gauge, the background gauge in the vacuum sector. 
C {x,y; j3) solve the transfer equations: 



among the coefficients with n > 1. Because 
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(4) 



The high-temperature expansions C^ {x, y; (3) = 
En=o'^tix,y)P'', 4ix,x) = I, trade the PDE Q by 
the recurrence relations 



[nl + {xk - yk){dkl - lQk)]ctix,y) - 
= [AI-Q±9,-y±]c±_i(f,y) 
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and Ci{x,x) — —V^{x), we obtain in the /3 = -infinite 
temperature- limit: 



ind2?=-^tr{ci(H+)-ci(7^-)} ^- [ d^xVi2i 
Att t^ J 
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the dimension of the self-dual vortex moduli space is 21. 
3. Standard lore in the semi-classical quantization of 
solitons tells us that the one-loop mass shift comes from 
the Casimir energy plus the contribution of the mass 
renormalization counter-terms: AMy = AMy -I- AMy. 
The vortex Casimir energy with respect to the vacuum 
Casimir energy is given formally by the formula: 



AMS = 



hm 



STr*{n+)' -STr(7^o)' 



where m = ew is the Higgs and vector boson mass at the 
critical point k = 1 . We choose a system of units where 
c = 1, but h has dimensions of length x mass. The 
"super traces" encode the ghost contribution to suppress 
the pure gauge oscillations: STr* {H+)^ = Tr* (H+)^ - 

Tr(7i:'^)^ and STr(Ho)^ = Tr(7i:o)^ - Tr(7i;^). The 
trace for the ghosts operators is purely functional: i.e., 
7i^ = - A +\'ip\'^, TLq = - a +1 are ordinary -non- 
matricial- Schrodinger operators. The star means that 
the 2n zero eigenvalues of 7i+ must be subtracted because 
zero modes only enter at two-loop order. 

In a minimal subtraction renormalization scheme, one 
adds the counter-terms £fj = /im/[|0p — l], £;?( = 



-iftm/yl^A'" with / 
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to cancel the di- 



vergences up to the one-loop-order that arises in the 
Higgs tadpole and two-point function, and in the two- 
point functions of the Goldstone and vector bosons. Fi- 
nite renormalizations are adjusted in such a way that 
the critical point k = 1 is reached at first-order in the 



loop expansion. Therefore, the contribution of the mass 
renormahzation counter-terms to the vortex mass is: 
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AM{; = AM^t^ 
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where S(V,T4) = / dx^ [(1 - l^p) - W^Vk]. 

We regularize both AMy and AMy by means of gen- 
eralized zeta functions. From the spectral resolution of 
a Fredholm operator 7i^„ = A„^„, one defines the gener- 
alized zeta function as the series C,n{s) = ^„ -^, which 
is a meromorphic function of the complex variable s. We 
can then hope that, despite their continuous spectra, our 
operators fits in this scheme, and write: 
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where Cwo(s) 



— ^^^ — nT) ^"^^ /i is a parameter 
of inverse length dimensions. Note that AMy — 
l[ia^^_iAM^{s), AMP = lim^^i AM^(s) and / == 



lim„ 



4 2m^L^ 



CWoW- 



4. Together with the high-temperature expansion the 
Mellin transform of the heat trace shows that 
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structure. Contrarily to the (l-l-l)-dimensional case, the 
value s = —^ for which we shall obtain the Casimir en- 
ergy is not a pole. 

Writing Cn = Yll=i[<^ri]aa{'H+) - Cn{W^), the contri- 
bution of the first coefficient to the Casimir energy 
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and exactly cancels the contribution of the mass renor- 
mahzation counter-terms -also finite for s = i-: 
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Subtracting the contribution of the 2/ zero modes we fi- 
nally obtain the following formula for the vortex mass 
shift: 
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is the sum of meromorphic and entire -Bu{s)- functions 
of s. Neglecting the entire parts and keeping a finite 
number of terms A^o in the asymptotic series for C'h(s), 
we find the following approximations for the generalized 
zeta functions concerning our problem: 
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7[s -I- n — 1, 1] — L dp (3''^" ^e '^ is the incomplete 
Gamma function, with a very well known meromorphic 



5. Computation of the coefficients of the asymptotic 
expansion is a difficult task; e.g. the second coefficient 

c+(f , ^) = -^ A V+{x) + l-Q+{x)Qt{x)V+{x) - 
-ldkQt{x)V+{x) + ^Q+{x)dkV+{x) + l[V+]\x) 

Defining the partial derivatives of the coefficients at y = 

X as 

y^x OXi 0X2 

we write their recurrence relations 
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starting from '^^''''>Cq^{x). 

We notice that [cn]jp{x) = (°'°)C^P(f) and thus 

Things are easier if we apply these formulae to cylindri- 
cally symmetric vortices. The ansatz 0(r, 9) — /(r)e''* 
and rAg{r,9) = la{r) plugged into the first-order equa- 
tions leads to: 
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Solutions of ^ with the boundary conditions lim f{r) — 

r — ^oo 

1, lim a{r) — 1, zeroes of the Higgs and vector fields at 

r — >oo 

the origin, /(O) = 0, a{0) — 0, and integer magnetic flux, 
^9 = — J ^ dOAg — 27r/, can be found by a mixture of 
analytical and numerical methods 13j . Henceforth, we 
shall focus on the case / = 1. 

The heat kernel coefficients depend on successive 
derivatives of the solution. This dependence can increase 
the error in the estimation of these coefficients because 
we handle an interpolating polynomial as the numerically 
generated solution, and the derivation of such a polyno- 
mial introduces inaccuracies. It is thus of crucial im- 
portance to use the first-order differential equations © 
in order to eliminate the derivatives of the solution and 
write the coefficients as expressions depending only on 
the fields. The recurrence formula now gives the coeffi- 
cients of the asymptotic expansion in terms of /(r) and 
air), e.g.: 
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E [^3]"(^' ^) = Y^^ [-4 "M' - 28 r^ a(rf (2 + 5 f{rf) + 

1 — 1 

+4:a(r)'^ ('20 + 9r2 + 'i2r'^ f(rf\ -2r'^ a(r) (-4 (l6-|-9r^) -f 
-I- (32 + 331 r^) f{rf + b7r^ f(r)'^\ +a{Tf (-256- 144 r^ 
-117r'' + 2r2 (56 + 183 r^) /(r)^ + 99r'' /(r)"*) + r* (-16+ 
+151 r^ + (392 -321 r^) f(rf + {-20 + l%9r'^) f{rf 
-29r^f{rf)] . 
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Plugging in these expressions the partially analytical 
partially numerical solution for /(r) and a{r), it is pos- 
sible to compute the coefficients -also for the ghost op- 
erator via similar but simpler formulae- and integrate 
numerically them in the whole plane. Thus, formula Q 
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provides us with the one-loop vortex mass shift, where 
we recall that 



J2[Cn]aain+) - Cnin^ 



The results are shown in the Table I: 



TABLE I: Seeley Coefficients and Mass Shift 
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EtiC(H+) 


c„(H^) 


iVo 


AMv{No)/hm 
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30.3513 


2.677510 


2 


-1.02814 
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13.0289 


0.270246 
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-1.08241 
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4.24732 


0.024586 
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-1.09191 
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1.05946 


0.001244 


5 


-1.09350 
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0.207369 


0.000013 


6 


-1.09373 



The final value for the vortex mass at one- loop order is: 

My = m f— - 1.09373ft) -I- o{h^). 

The convergence up to the sixth order in the asymptotic 
expansion is very good. We have no means, however, of 
estimating the error. In the case of A((?!))| kinks we found 
agreement between the result obtained by this method 
and the exact result up to the fourth decimal figure, see 
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